Spatial patterns in optical lattices submitted to gauge potentials 
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Abstract. - We study the vortex formation in optical lattices submitted to artificial gauge po- 
tentials. We compute the superfluid density for Abelian and non-Abelian gauge potentials with a 
mean-field approach of the Bose-Hubbard model and we determine the rule describing the num- 
ber of vortices as a function of the effective magnetic flux. This simple rule is represented by 
a remarkably rich figure that represents the superfluid density as a function of the flux. The 
phenomena which emanate from this work should be observed experimentally in optical lattices 
within which atom tunneling is laser-assisted and described by commutative or non- commutative 
tunneling operators. 



The observation and comprehension of spatial patterns 
in physical systems constitute a major field of investigation 
[1]. The study of these complex structures, which arise at 
the macroscopic level, is necessary in order to understand 
the behaviour of numerous quantum systems such as type- 
II superconductors [2-5], rotating helium superfiuids [6,7], 
rotating Bose-Einstein condensates [8,9] and the fractional 
quantum Hall effect [10]. 

In the context of superfluid Bose gases [8,9], these ob- 
servable patterns result from topological defects in the 
macroscopic wave function which constitute an order pa- 
rameter for the insulator-superfluid transition [11, 12]. 
These singular points in which the superfluid density dra- 
matically decreases, known as vortices, are regarded as 
interacting entities and usually form regular lattices. 

The very recent progresses in cooling and manipula- 
tion of cold atomic systems [13-16] have open the way to 
investigations of fundamental aspects of condensed mat- 
ter physics. Optical lattices, created by adjusted lasers, 
constitute an elegant way to study periodically trapped 
atomic systems [17,18]. Among the various experimental 
configurations that can now be set up is the application 
of artificial gauge potentials (artificial magnetic fields) to 
neutral atoms [17,20-23]. Such setups notably offer the 
possibility to reproduce the dynamics of periodically con- 
strained electrons submitted to a magnetic field [17,20]. 
This should lead to the observation of the corresponding 
fractal energy spectrum (the Hofstadter butterfly [19]) and 
to the exploration of the quantum Hall effects [21,24]. In 



particular the realization of non-Abelian gauge potentials 
has been envisaged [17], allowing the observation of a non- 
Abelian Aharonov-Bohm effect [22], magnetic monopoles 
[23], particular metal-insulator transitions [25] and an in- 
teger quantum Hall effect for neutral currents [24]. 

In this context, a problem of interest is the effect of 
gauge potentials on the Mott insulator transition: It is 
known that the insulator-superfluid transition is affected 
by an Abelian gauge potential [26]: the Mott lobes, which 
characterize the corresponding phase diagram, are in- 
creased when the artificial magnetic field is turned on. 
Moreover the critical point of the transition is proportional 
to the bandwidth of the Hofstadter butterfly [19], the en- 
ergy spectrum of the noninteracting problem. 

In this work we are interested by the effects induced by 
a gauge potential once the system has reached the super- 
fluid regime. Since the creation of an effective magnetic 
field within an optical lattice reproduces the dynamics of 
a spinning Bose system, as expected we observed in this 
work vortex configurations [8,9,26-28] that we studied for 
general non-Abelian gauge potentials. In this letter we 
first introduce the Bose-Hubbard Hamiltonian describing 
the optical lattice submitted to non-Abelian gauge po- 
tentials and we briefly describe the mean-field approach 
considered in this context. The vortex configurations are 
studied for Abelian and non-Abelian configurations of the 
external gauge potential. The general rule describing the 
pattern formation that emanates from this work is even- 
tually exposed. 
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To start with, we consider a bosonic system trapped in 
a two-dimensional optical lattice and submitted to artifi- 
cial gauge potentials. In this letter we consider a family of 
non-Abelian gauge potentials, which contains Abelian po- 
tentials as particular cases. The system is then described 
by a generalized Bose-Hubbard Hamiltonian written as a 
sum of single-site terms H = ^ 'Hn,m^ with 

Wn,m = —ta{cili,m\^x\ciri-\-l,m} — tb{cili^m\Uy\cLn,m-\-l} 
— ta{cili^fyi\U^\(ln—l,m } - tb{ai^jUl\an,m-i} 

+ -^({4,™!"",™} - l){4,ml«n,m} " ^{4,™!"",™} 

(1) 

The first two lines of the expression determine the dy- 
namics of noninteracting "charged bosons" in a two- 
dimensional lattice within a tight-binding approximation 
as studied experimentally in cold atom systems trapped in 
optical lattices [13-16]. For Ux^y belonging to the Abelian 
group of unitary complex numbers U{1), this noninteract- 
ing part of the Hamiltonian reproduces the evolution of 
an electronic system constrained by a periodic potential 
and submitted to a magnetic field [20]. In the present 
context one consider a system with U (2) gauge structure, 
which implies that the operators Ux and Uy are 2x2 uni- 
tary matrices and the single-particle wave function has 
two components: when tunneling occurs along the x (y) 
direction, the matrix Ux (Uy) acts as a tunneling operator 
on the two- component wave function and the system is 
non-Abelian for [Ux^Uy] ^ [22,25]. The operators Ux,y 
are related to the gauge potential A present in the system 
via the usual relation Ux,y — e*"^^'^. In order to indicate 
multicomponent wave functions and operators, we use the 
following notations introduced in Refs. [24,29]: A row of 
2 orthogonal kets is denoted by Wi))\ = (|'^)i, |'0)2) and 
a column of 2 kets by ||'0)}.The coefficients ta and tb de- 
scribe transfers in the x = nix and y = mly directions, 
where {Ix^ ly) are the unit cell's lengths, and ^} (resp. 
\(^n,m}) are the 2-component bosonic creation (resp. an- 
nihilation) field operators on site {n^m). The second line 
of the expression includes the interaction term and the 
chemical potential. 

The system is solved using a mean-field approach based 
on the decoupling formula [11,26,27] 

- (4,m,i)K,m,i') (2) 

with j,/ = 1,2. In the Abelian case, the average value 
{(^n,m,i) = ^n,m,i = ^n,m,2 defines an order parameter for 
the insulator-superfiuid transition and is non-zero when 
the ground state is in the superfiuid state [26]. As we are 
dealing with a non-Abelian system for which ^n,m,i gener- 
ally differs from ^n,m,2 we define |^n,mP = {^n,m|^n,m} 

as the superfiuid density, and refer |^n,m,j 

r as a "color" 

density with j = 1,2. 



We and others have recently studied the effects in- 
duced by non-Abelian gauge potentials in optical lattices 
[22,24,25]. Each vertex of a 3D optical lattice contains an 
atom with degenerate Zeeman sublevels in the hyperfine 
ground state manifolds {\g)j^ with j = 1,2, so that 
the states we are dealing with are represented by the rows 
{|e)| = (|e)i,|e)2) or {|^)| = 1^)2). Lasers are used 

in order to create the lattice, as well as the non-Abelian 
gauge potential through state-dependent control of hop- 
pings which take place within every plane z = constant. 
In such setups, atoms hop from a site to another and their 
tunnelings are described by unitary non-commutative op- 
erators [22]. The corresponding gauge potential is given 

by 

/ 2'imai 
I 27rna2 , 



A = 



,0 



(3) 



and induces an effective "magnetic" field characterized by 
two effective magnetic fluxes ai and ce2, which can be 
controlled with external lasers [22]. This gauge poten- 
tial is Abelian under the condition that o^i — a2 is an 
integer and non-Abelian otherwise. The gauge poten- 
tial preserves the translational symmetry along the 7/-axis, 
such that |^n,m} = l^n}- We consider quantized fluxes 
= Pi/q and 0^2 = withpi,p2 5 Q integers, such that 
l^n+q} = l^n}- Including the potential A in the mean- 
field Hamiltonian (1,2), one finds the non-zero matrix el- 
ements of in the occupation number basis \Nnj)^ 
with Nnj = 0, 00: 



V 



{Nn,l\K'' \Nn,l) = -(A^n,l " l)A^n,l " /iA^n,l 
+ ta(^;,l^n+l,2 + ^;,2^n+l,l + h.C.) 

-ht6(2cos(27rQ;in)^; + 2 C05(27ra2ri)^;,2^n,2) 



{Nn,l + l|Wr |^n,l) = -taVA^n,l + l(^n+l,2 + ^n-1,2) 

- 2t6V']Vn7TTco5(27rain)^n,i 



(A^n,2 + l\n^^\Nn,2) = "^a V^n,2 + 1 (^n+1,1 + ^n-l,l) 



2tb^/N^^iJTlcos{27^a2n)^r^,2 



(4) 



where the expressions anj\Nnj + 1) = -s/Nnj + l|A^n,i) 
and al^ j\Nnj) = \/Nnj + M^nj + 1) have been used. 
The ground state \Gn) is computed in a truncated basis, 
Nnj = 1, Nmax, and the order parameter 1^^} is given 
by = {Gn\ctn,j\Gn)- The latter expressions form a set 
of self-consistent equations which is solved in order to find 
the ground state \Gn) and the order parameter |^n}- 

The system presents symmetries under translations 
along the ?/-axis and the system is g-periodic along the 
X-axis when the gauge potential (3) is turned on: one com- 
pute the densities on a supercell qxl with periodic bound- 
ary conditions. The superfiuid density and the "color" 
densities are depicted in Fig.l for specific Abelian (Fig.l.a) 
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Fig. 1: "Color" densities |^n,i|^ (red line) and |^n,2|^ (black 
line) and superfluid density |^n|^ = {^n|^n} (thick blue line) 
for /i 1.4, ta = h = 0.05 and q 95: (a) ai = a2 = |; (b) 
tti = ^ and ^2 = I 

and non-Abelian (Fig.l.b) situations when the system is in 
the superfluid regime (/i = 1.4, ta = h = 0.05). For both 
cases oscillations in the superfluid densities are observed; 
in the following we consider the local minima as vortices. 
In Fig.l, one set the supercell length to = 95 and one 
finds 3 vortices when = 0^2 = | (Fig.l. a). In the non- 
Abelian case (Fig.l.b) one finds 7 vortices when ai = ^ 
and a2 = I . We notice that the number of vortices relates 
to the values of pi and p2 and in the non-Abelian case we 
notice that the value of pj mainly influences |^n,jP in the 
sense that |^n,iP (red line) has 7 vortices while |^n,2p 
(black line) has 3 main vortices. We also notice that the 
general shape of the superfluid density (thick blue line) 
shows 3 principal slopes. 

A simple and general rule is observed for arbitrary val- 
ues of the fluxes ai and a2. In the Abelian case, for 
which ai = 0L2 = p/q^ with p < one naturally finds 
that |^n,iP = |^n,2p and exactly p vortices are formed 



in the supercell of length q (see Fig.l. a). In the non- 
Abelian case, for which pi 7^ p2, the "color" densities differ 
|^n,iP 7^ |^n,2p and the situation is complex: one finds 
p vortices, with p = max{pi^p2)^ but the overall shape 
of the total density |^nP is mostly influenced by with 
p = min{pi^p2)'. in Fig.l.b one can seep = 7 local vortices 
and p = 3 general slopes. 

This important result leads to a striking representation 
of the superfluid density as a fonction of the fluxes. In 
Fig. 2. a which illustrates the Abelian case, one observes 
the manner by which the number of vortices (yellow dots) 
evolves with the flux a = ai = a2'- there are exactly p 
vortices within the supercell at a = |, for p < q/2^ and 
the superfluid density's general shape smoothly follows the 
rule through concentric circles (see the contour plot in 
Fig. 2. a). We notice that the figure is symmetric around 
the value a = 0.5 and is periodic in a with period one. 
In order to represent the non-Abelian case (Fig.2.b) one 
set Q^i = ^ ~ and q = 149. In this case, ai ~ ^ 
with pi = 15 and we note that the result is robust for 
ai irrational. The contour plot shows that the superfluid 
density is mainly affected by the fixed value pi = 15 which 
generates 15 general slopes in the density for every value of 
P2 > 15, but a pattern similar to the analoguous Abelian 
case is still observable in the background which correspond 
to the many p = max{pi,p2) vortices formed as decribed 
here above. 

In the Abelian case, the superfluid density's represen- 
tation in terms of the flux shows a particular property 
which we refer to as "scale invariance". This property is 
summarized as follows: when the flux is rational, a = -, 

' q ' 

we have shown that p vortices are formed whithin a (7 x 1 
supercell; but when considering a sample of Q sites within 
the supercell, one finds exactly p vortices at a ~ This 
property is illustrated in Fig. 3 which shows the superfluid 
density computed in a 109 x 1 supercell and represented 
as a function of the flux in a sample of length Q = 17: 
as expected p vortices are formed at values a = j^. An 
important consequence of this result is the fact that the 
number of vortices is independant of the supercell length's 
value q. Moreover, this property gives rise to the complex 
structure represented in Fig. 2a: the structure is formed 
in such a way that, between any interval n = [1,(5] of 
arbitrary length Q, exactly p dots are drawn at a = ^. 

The previous results presented here above are obtained 
for specific values of the hopping parameters ta^ h and the 
chemical potential fi. In Fig. 4. a we show the superfluid- 
insulator phase diagram in the jj. — t plane for the non- 
Abelian case, with t = ta — thi Pi = 1, P2 = 4 and q = 48: 
three distinctive Mott lobes, inside which l^nP = 0, 
are clearly observed. In Fig.4.b we represent the super- 
fluid density while progressing vertically in the insulator- 
superfluid phase diagram at t = 0.05: the superfluid 
density |^nP is vanishing in the whole supercell around 
/i = 0.5 which corresponds to the entrance of the first Mott 
lobe, then in the superfluid region, the vortex formation 
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Fig. 2: Superfluid density |^n|^ = {^n|^n} as a fonction of a2 
(2/-axis) for /i = 1.4, ta = tb = 0.05: (a) ai — a2 with q — 109 
(yellow dots represent the vortices); (b) ai = ^ ~ with 
q — 149 and n = 1, ^ (x-axis) 



is important when approaching the Mott lobes at /i = 1.5 
and /i = 2.5. The same phenomenon is observed in the 
Abelian case illustrated in Fig. 5 for a = ai =0^2 = 
four vortices are formed around /i = 1.5 (near the second 
Mott lobe) and /i = 2.5 (near the third Mott lobe) and 
vanish around /i = 0.5 (first Mott lobe). 

In summary, the pattern formation resulting from gauge 
potentials created in optical lattices have been studied. 
The number of vortices formed in the superfluid Bose 
gase directly relates to the gauge potential values. In 
the Abelian case, this number results from a simple rule 
which leads to a salient figure that illustrates the super- 
fluid density as a function of the effective magnetic flux. 
This phenomenon is observed in a wide region of the su- 
perfluid regime. As the vortices studied in this context 
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Fig. 3: Superfluid density |^n|^ = {^n|^n} as a fonction of 
a2 (y-axis) for /j, = lA, ta = tb = 0.05, ai = a2 — ^ for a 
supercell of length q = 109 and n = 1, .., 18 (x-axis). Yellow 
dots represent the vortices and the thin blue lines show the 
constant functions = = 77, and ^ where 1, 2 and 3 
vortices are respectively observed between n = 1 and n — Q — 
17 

have typical size L ^ b — 1{) (with typical lattice con- 
stant ao = 0.5 — 5 /im, see Refs. [12,30]), the patterns 
discussed in this work should be directly observed in lab- 
oratories by optical methods as suggested by Wu et al. in 
Ref. [27]. 
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Fig. 5: Superfluid density |^n|^ = {^n|^n} as a fonction of the 
chemical potential fi (^/-axis) for ta — th — 0.05, oli — a2 — ^, 
q — AS and n — 1, .., 48 (x-axis) 
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